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Abstract 

We use the mesoscopic nonequilibrium thermodynamics theory to de- 
rive the general kinetic equation of a system in the presence of potential 
barriers. The result is applied to the description of the evolution of sys- 
tems whose dynamics is influenced by entropic barriers. We analyze in 
detail the case of diffusion in a domain of irregular geometry in which the 
presence of the boundaries induces an entropy barrier when approaching 
the exact dynamics by a coarsening of the description. The corresponding 
kinetic equation, named Pick- Jacobs equation, is obtained, and its valid- 
ity is generalized through the formulation of a scaling law for the diffusion 
coefficient which depends on the shape of the boundaries. The method we 
propose can be useful to analyze the dynamics of systems at the nanoscale 
where the presence of entropy barriers is a common feature. 

1 Introduction 

The free energy landscape of a complex system, when represented as a func- 
tion of an order parameter or reaction coordinate, presents an intricate aspect 
consisting of multiple local minima, defining metastable states, separated by 
barriers. The nature of the barriers depends on which thermodynamic poten- 
tial varies when passing from one well to the other and their presence plays an 
important role in the dynamics of the system. Whereas energy barriers are more 
frequent in problems of solid state physics, entropy barriers are often encoun- 
tered in soft condensed matter and biological systems. These barriers may also 
appear when coarsening the description of a complex system in order to simplify 
its dynamics. The elimination of some coordinates can be performed by intro- 
ducing a configurational entropy accounting for the degeneracy (i.e. the number 
of equivalent configurations) of the state defined by the remaining coordinates. 
This procedure has been used for instance in the context of glasses 0]. 

The study of the dynamic properties of the system at the mesoscopic level 
requires the knowledge of the probability distribution function of the relevant 
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degrees of freedom which evolves according to a kinetic equation of the Fokker- 
Planck type. This equation has been usually derived from kinetic theory in the 
diffusion limit |2|-|Q|, from a master equation jH), and from projection operator 
techniques |§, |^ | ■ 

Our purpose in this paper is to propose a new and simple theoretical frame- 
work to analyze the mesoscopic dynamics of systems in the presence of energetic 
and entropic barriers. The procedure we follow is based on the application of 
mesoscopic nonequilibrium thermodynamics and has been previously 

used to analyze transport and activated processes in systems governed by a dy- 
namics of the Fokker-Planck type. The kinetic equation follows after obtaining 
the probability current, occurring in the corresponding continuity equation in 
the configuration space, from the entropy production accounting from dissipa- 
tion in that space related to the underlying diffusion process of the probability 
density. That equation applies irrespectively of the nature of the equihbrium 
state of the system. When the system is in contact with a heat bath and its 
volume remains unaltered, in which case the proper thermodynamic potential is 
the free energy, our equation describes the dynamics in the presence of energetic 
and entropic barriers. The validity of our method, however, goes beyond this 
standard situation and embraces other cases of interest characterized by other 
thermodynamic potentials and statistical ensembles. 

The paper is distributed in the following way. In Section II, we present 
a general derivation of the kinetic equation in the framework of mesoscopic 
nonequilibrium thermodynamics. That equation governs the evolution of the 
probability density in the presence of barriers of any nature. In Section III we 
discuss an example in which the presence of an entropic barrier manifests: dif- 
fusion in a channel of varying cross section. The corresponding kinetic equation 
has been referred to as the Fick- Jacobs equation. Section IV is devoted to the 
derivation of that equation from mesoscopic nonequilibrium thermodynamics 
and to the formulation of a scaling law for the diffusion coefficient which is 
compared to previous results. In Section V we extend the results of Section IV 
to the case in which an external field acts on the diffusing particle. Finally, in 
Section VI we present our main conclusions. 



2 Kinetic equations in the presence of potential 
barriers 

Our purpose in this section is to present a general derivation of the kinetic 
equations describing the evolution in time of the probability density when the 
statics of the system is characterized by a thermodynamic potential landscape. 

We consider that the equihbrium state of the system is described by the 
probability distribution function 

Pe,-exp(^^^j, (1) 

where W(x) is the minimum reversible work required to change the state of 
the system which is related to the maximum useful work which can be 
extracted from it, sometimes referred to as available energy ||2^ or availability 
k is Boltzmann's constant, T is the temperature of the medium and x an 
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arbitrary set of coordinates which may represent the velocity of a particle, the 
orientation of a spin, the size of a macromolecule or whatever coordinate or 
order parameter whose values defines the state of the system in a phase space. 
Variations of the minimum work for a thermodynamic system are expressed as 

AW = AC/ - TAS + pAV - ^lAN + yAY + . . . , (2) 

where U is the internal energy, S the entropy, V the volume, and N the number 
of particles of the system, whereas T, p, and fi are the temperature, pressure and 
chemical potential of the environment. The term yAY represent other kinds of 
work (electric, magnetic, surface work...) performed on the system, being y the 
intensive parameter and Y its conjugate extensive variable ||2^. For the case 
of constant temperature, volume and number of particles, the minimum work 
corresponds to the Helmholtz free energy F. In general, that quantity reduces to 
the thermodynamics potentials by imposing to expression (||) the corresponding 
constraints ||2^ . 

Having specified the statics of the system we will derive the kinetic equation 
describing the evolution of the probability density. To this purpose, we will 
use the framework of mesoscopic nonequilibrium thermodynamics. This theory 
applies the scheme of nonequilibrium thermodynamics to describe the dynamics 
of mesoscopic degrees of freedom. The treatment of a diffusion process in the 
framework of nonequilibrium thermodynamics is extended to the case in which 
the relevant quantity is a probability density, defined in phase space, instead of 
a mass density. The starting point is then the formulation of the Gibbs equation 

6S=-^ J fi{x)SP{x,t)dx, (3) 

which resembles the corresponding law proposed in nonequilibrium thermody- 
namics in terms of the mass density of particles. Here fi{x) is a generalized 
chemical potential whereas P{x, t) is the probability density. 
This expression is compatible with the Gibbs entropy 

S{t) ^-k f Fix, t) In ^^^^dx + Seg , (4) 

where Seq is the entropy when the system and the heat bath are at equilibrium. 
Effectively, taking variations in Eq. (^) and using the expression of the chemical 
potential per particle |^ 

= kT In + ^leq , (5) 

where ^eq is the chemical potential at equilibrium, one arrives at expression @. 

From the Gibbs equation we can obtain the entropy production related to 
the underlying diffusion process in x-space. Calculating the time derivative of 
the entropy from Eq. one has 

dS , f , F dF , 

Quite generally, we may assume that the probability density satisfies the conti- 
nuity equation 
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where J(x, t) is a current defined in x-space which has to be determined. Sub- 
stituting now the time derivative of Eq. into Eq. (||) one obtains the entropy 
production 



f d P 
'k I Jix, t) ■ — In dx . 

J ox Peg 



(8) 



To arrive at Eq. (g) we have also performed a partial integration and assumed 
that the current vanishes at the boundaries of the system in x-space. 

In the nonequilibrium thermodynamics scheme, the entropy production 
consists of contributions of products between the current J_{x, t) and its con- 
jugated thermodynamic force — fc^ln-^, for each value of the coordinate x. 
Assuming locality in x-space, for which only currents and forces at the same 
value of x are coupled, one obtains the linear law 

^=-fcLfe)|ln^, (9) 

in which the phenomenological coefiicient L{x) may in general depend on x_. 

The resulting kinetic equation then follows by substituting Eq. (||) into the 
continuity equation (|^) 

dP d f d , P 



at dx v^^a^^^^^)' ^''^ 



where we have defined the diffusion coefficient 

^fe) = ^. (11) 
This equation, which in view of Eq. ([l]) can also be written as 
dP_d_ f dP D_dAW \ 

dt ^ dx' [ dx + kT dx ) ' ^ ' 

is the Fokker-Planck equation accounting for the evolution of the probability 
density in a;-space. 

Under the conditions for which W = F = U — TS, this equation transforms 
into the Fokker-Planck equation for a system in the presence of energy and 
entropy barriers. One then obtains 

at dx \ dx kT dx k dx J ' ^ ' 

where the drift term consists of contributions due to an external potential and 
to variations of the entropy. When the barrier is purely energetic, only the 
first contribution remains, and the previous equation reduces to the well-known 
Fokker-Planck equation 

dt ^ dx' [ dx ^ kT dx r ^ ' 
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If the nature of the barrier is purely entropic, Eq. (13) then corresponds to the 
Fick-Jacobs equation 



dP 
'dt 



d_ 

dx 



( 




k dx 



P 



) 



(15) 



first derived in the context of diffusion of Brownian particles in a channel of 
non-constant cross section P3|. 

The general form of Eq. (|l2|) , in which the equilibrium distribution function 
does not need to be specified and is given in general by Eq. (|l|), makes that 
result appHcable to a great diversity of situations. For example, for a Brownian 
particle, for which the minimum work is simply its kinetic energy, that equation 
correspond to the usual Fokker-Planck equation. When x_ is an order parameter 
or a reaction coordinate, Eq. (|2|) provides the kinetic equation in the presence 
of barriers. The method used in this section then offers a common formalism 
able to analyze the dynamics of a system in the presence of energy and entropy 
barriers. In this paper, we will mainly focus our analysis to the case in which 
the potential is strictly entropic. 

3 Diffusion through a channel of varying cross 
section: the Fick-Jacobs equation 

The infiuence that the presence of an entropic barrier exerts on the dynamics of a 
system can be illustrated by deriving the kinetic equation governing Brownian 
diffusion through a channel of varying cross-section. This equation was first 
proposed by .lacobs ||2^ and subsequently rederived by Zwanzig ||2^ on the 
basis of more fundamental arguments. 

3.1 Jacobs' derivation 

Jacobs in his book Diffusion Processes ||2^ provided an heuristic derivation of 
the equation governing diffusion in a symmetric tube whose cross section A{x) 
varies along the x axis, defined by the center line of the tube. The argument 
runs as follows. 

Consider an elementary volume of thickness dx perpendicular to the axis of 
the tube. The total amount of particles in this slice, at x and time t is C{x, t)dx, 
which is the integral of the concentration over the volume A{x)dx. The rates 
of entrance and exit of the diffusing substance into this volume are given by 
the Fick's law J = — Do^§^, where Dq is the diffusion coefficient, and C/A 
is the local volume concentration. In this case, both rates are different not 
only because the concentration gradient |^ changes with the distance, but also 
due to the variation of the cross-section of the channel. Explicitly, the rate of 
entrance into this elementary volume is —DqA ^^^J^ , whereas the rate of exit 



is given by —Do yA—^^ — h ^ \^A j dx + ...j . The difference between both 
rates provides the rate of change of the substance in the elementary volume, 
which can also be expressed as ^dx. Neglecting quadratic terms in dx, one 
easily arrives at the equation governing the diffusion in the channel 
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9^ 9 f . d fC\\ ^ d fdC 1 dA{x)\ 

-m - ^'Yx r 9^ \a))- ^'Yx [j^ - Mx)^^ J ^ ^ 

which is referred to as the Fick- Jacobs equation. 
3.2 Zwanzig's derivation 

Zwanzig reported in Ref. a more general and rigorous derivation of the 

Fick-Jacobs equation. For the sake of simplicity, let us review the 2D case. The 
method consists of performing a reduction in the number of coordinates from 
the 2D Smoluchowski equation to a ID description. 

The starting point is the two-dimensional Smoluchowski equation for diffu- 
sion through a general potential U{x, y) 

^A^gA = Do-^e-^^^^-y^^e^^^--y^cix,y,t) (17) 
+Z?o^e-'5^("'?')^e'^^("'«)c(x,2/,i) (18) 

where P — 1/kT, and c(x,y,t) is the concentration. To perform the reduction 
to ID, this equation is integrated over the variable y, leading to 

dC{x,t) a I ^^m^,y)d,^i.,y)^(^^^y^t)ay (19) 



dt dx J dx 

where the reduced concentration is defined as 



Cix,t) = j cix,y,t)dy. (20) 

The key point of the derivation is the assumption of equilibration in the 
transverse direction. Under this assumption, one can define an averaged x- 
dependent free energy F{x) through the expression 

e-/3F(x) ^ fe-m.-.y)dy (21) 



from which one can define a normalized conditional probability distribution 

^~PU{x,y) 

Then, under the local equilibrium approximation, the concentration c{x,y,t) 
factorizes as follows 

c{x,y,t)^C{x,t)p{x;y) (23) 
Taking these considerations into account, one finally obtains 

^ ^D^e-P^(^^ |_e/3^^(-)C(x, t) (24) 

which constitutes the generalization of the Fick-Jacobs equation ( [T^ ) to the case 
of a two dimensional potential U{x,y). 
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In the previous analysis, we have not taken into account the fact that the 
normal flux must vanish at the boundaries. The role played by this zero normal 
flux conditions, can be replaced by the conflning potential U{x,y). In fact, if 
U{x, y) is a box-like potential, i.e. zero inside the tube and infinite outside, it is 
clear that e~^^'--^^ = 2w{x), being w{x) the half width of the tube. In this case, 
the barrier is purely entropic, and one recovers the usual Fick- Jacobs equation. 

The extension to the three-dimensional case is quite straightforward by tak- 
ing into account that in that case the integration in the transverse coordinates 
involves variables y and z, and the width of the 2D tube must be replaced by 
the transverse section nw{x)'^ . 

Zwanzig's analysis clearly manifests that the accuracy of the Fick-Jacobs 
equation is conditioned to the existence of local equilibrium in the transverse 
direction. The author analyzed the effect of the deviations from the local equilib- 
rium related to the variations of the density Sc{x, y, t) = c(x, y, t)—C{x,t)p{x; y). 
He derived an equation for the evolution of this deviations, which suggested that 
the accuracy of the Fick-.Iacobs equations is restricted to situations verifying 
|w'(a;)| < 1, that is, when the section of the tube varies smoothly. In addition, 
he showed that the range of validity of the one-dimensional Fick-Jacobs descrip- 
tion could be extended by introducing a position-dependent effective diffusion 
coefficient 

1 

' 1 + ^w'{xY 

where Dq is the molecular diffusion coefficient and the parameter 7 depends 
on the dimensionality. The explicit expression for this coefficient is obtained 
through an expansion in powers of w'{x) iQ. The result is 

D{x)^DQ{l--iw'{xf + ...) (26) 

where 

being v the transverse coordinate scaled by the function w{x) (i.e. v = 
in 2D, and v — for a 3D tube of radius r with cyHndrical symmetry), 

and V{v) — U{x,y). For a purely confining potential, the results obtained by 
Zwanzig are 7 = 1/3 for the 2D case, and 7 = 1/2 for a 3D tube with cylindrical 



Dzix)^Do ^ , (25) 



symmetry. From the expansion ( |26| ) he infers the expression (|25| ) based merely 
on the fact that this reconstruction of the series improves the agreement with 
the exact results. 



4 The Fick-Jacobs equation from mesoscopic nonequi- 
librium thermodynamics 

For the case of an enclosure of varying cross-section, the concept of entropic 
barrier is remarkably simple. At equilibrium, the density of diffusing material 
po is constant. If we contract the 3D description retaining the coordinate x, the 
resulting ID equilibrium distribution is 
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Peq{x) = J podydz ^ poA{x) . (28) 

In this case, the diffusing particle constitute an isolated system. Therefore 
its corresponding thermodynamic potential obtained form Eq. (||) is simply 
AW = —TAS, where the entropic barrier is, in accordance with Eq. (||), 

AS{x) ^ kliiA{x) . (29) 

The previous equation clearly manifests that the entropic barrier originates 
from the variation of the space available for the diffusing particles. Notice that 
Eq. ( pol) corresponds to the usual microcanonical definition of the entropy in 
terms of the number of states which in this case is simply proportional to the 
area of the tube. 

Once identified the equihbrium distribution (or the entropic barrier) , the ki- 
netic equation describing diffusion in the presence of the entropic barrier follows 
from our general scheme developed in Section 2. According to Eq. (|l5|), the 
kinetic equation is 

which has the same structure as the Fick- Jacobs equation, but with a spatial 
dependent diffusion coefficient. 

At this point, it is worth to analyze the role played by this coefficient. In Ref. 
pi] it was shown that the validity of the ID description could be extended if we 
use an effective spatial-dependent diffusion coefficient. However, the agreement 
one can achieve with the expression (|5|) proposed by Zwanzig is not always 
satisfactory. In particular, for the 3D example discussed in Ref. ||2^, results 
for the current of particles obtained by using the expression of the diffusion 
coefficient ( p5|) present severe discrepancies with the corresponding expression 
coming from the exact result of the 3D Smoluchowski equation, when the section 
of the tube changes abruptly. 

The appearance of an entropic barrier originates from the reduction of the 
space to a single coordinate. But this reduction may also have implications on 
the form of the diffusion coefficient. The molecular diffusion coefficient in the 
real space, Dq, gives information about the dispersion of the displacement. In 
2D, we can estimate 

which manifests that the diffusion coefficient involves the displacement in both 
coordinates. Once we have contracted the y coordinate, the resulting effective 
diffusion coefficient only depends on the dispersion in the remaining x coordi- 
nate, that is qualitatively Deff ^'^^^ . Therefore, from Eq. (^l|) we can infer 
the behaviour D^j / ~ ^i^l^-^y ' Previous equation provides a hint about the 

fact that the reduction of coordinates may involve not only the appearance of 
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an entropic barrier, but also a scaling of the diffusion coefficient. Following this 
heuristic reasoning, we will propose the scaling law 



^^(^) = ^°(TT7(^ (32) 

for the diffusion coefficient appearing in Eq. (|30|), where y'{x) = and y{x) = 
w{x) defines the shape of the enclosure. The objective is then to test if the use of 
this expression is able to extend the range of validity of the ID Fick-Jacob-type 
description. 

The values of the scaling exponents can be estimated from the calculations 
performed in Ref. iQ. We can expand our expression for the diffusion coefficient 
Ds{x) in terms of and compare the result with the one obtained by 

Zwanzig (Eq. From that comparison, one obtains that a reasonable choice 
of the scaling exponents is a = 1/3 for the 2D case and a — 1/2 for the case of 
a 3D tube. 

In the remaining of this section, we will test this result for the two cases 
discussed in Ref. |Q which admit an exact solution. We will see that the 
agreement with the exact result, when we use our effective diffusion coefficient 
given by Eq. (^2|), is significantly improved. 

4.1 Effective Diffusion CoefHcient in a 2D periodic sym- 
metric channel 

The first case under scrutiny is the diffusion in a 2D periodic channel defined 

by 

X = M + a coshu sinu; y = v + asmhv cosu (33) 

where diffusion occurs in the region comprised between < w < 27r and —V< 
V <V. The parameter A = a cosh y must be smaller than 1, to avoid double 
valuation of the walls. Fig. |l| illustrates the shape of the tube for V = 0.5 and 
a = 0.5. 

At very long times, diffusion in periodic channels can be described through 
an effective diffusion coefficient D* , accounting for the effects of passing through 
many constrictions in the channel. This effective coefficient can be evaluated as 

il 

where ( ) denotes average over one period of the potential. 

For this particular geometry, an exact expression for D* was derived in Ref. 

11 

'D* 



Do 



1 



A^tanhy ^ 



V 



(35) 



and was compared to the results obtained with the Fick- Jacobs equation using 
a constant diffusion coefficient Dq, and with the position-dependent diffusion 
coefficient Dz{x). 

In Fig. 1^ we have represented the results for D* obtained using the different 
approaches mentioned above, plus the ones corresponding to our expression 
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Figure 1: Shape of the two-dimensional periodic channel defined parametrically 
by Eq. (H), for V = 0.5 and a = 0.5. 
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Figure 2: Effective diffusion coefficient D* for the 2D channel for V — 0.5 as a 
function of the parameter a, using different expressions for the diffusion coeffi- 
cient. The inset is a zoom to illustrate the accuracy of the different approaches. 



Ds{x). One can easily realize that the results obtained with the expression 
of Ds{x) we propose improve the ones when using Dz{x) and constitute an 
excellent approximation of the exact result. 

4.2 Steady-state flux through a 3D hyperboloidal cone 

The improvement of the results one can achieve with Dg (x) is more remarkable 
for the 3D example discussed in Ref. ||24). In this case, the shape of the tube 
corresponds to a 3D hyperboloidal cone as the one depicted in Fig. ^. It can be 
conveniently described by using oblate spheroidal coordinates (^,77, c/?), related 
to the cyHndrical coordinates (r, z,ip) by 

r'^ =a^{^^ + z^a^Tj (36) 

In this coordinate system, ^ = corresponds to the small hole and ^ = cxd to 
the far end of the tube; 77 = 1 is the z-axis, 77 = is the (x,y)-plane, whereas 
rj = r]o corresponds to the surface of the tube. Diffusion then takes place in the 
region 
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Figure 3: 3D hyperboloidal cone for a=l,0<^<5 and 770 = 0.5. 




Figure 4: Comparison of steady state fluxes through an arbitrary hyperboloidal 
cone, against the value of 770. 



< ^ < 00, ?]o <ri < 1, < < 27r (37) 

The steady-state 3D diffusion equation can be solved by using the boundary 
conditions: C = for ^ = 0, C = Co for ^ = 00, and the condition that the 
normal flux vanishes at the walls = rip. In this situation, the exact steady-state 
flux through the exit hole ^ = is 

Jexact = 4Da(l - ?7o)Cn (38) 

Alternatively, one can use the Pick- Jacobs equation derived in the previous 
section to calculate that flux, yielding 

''^'{r'^immr '''' 

The results are depicted in Fig. ^, corresponding to the "exact", the unmodifled 
Fick-Jacobs equation, in which D{z) = Dq is a constant, and to the choices 
Dz{z) and Dsiz). 

Prom the examples we have discussed, it becomes clear that our proposal 
of a scaling law for the spatial-dependent diffusion coefflcient leads to a better 
agreement with the exact results. 
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Figure 5: Schematic representation of the 2D tube in the presence of gravity. 



5 Diffusion in a channel of varying cross-section 
under the influence of a field 

In the previous section, we have discussed the effect of entropic barriers in the 
dynamics. The presence of that barriers is sometimes accompanied by energy 
barriers. The theory we have developed admits a generaUzation to the case 
in which both an entropic and an energetic barrier coexist. This will be the 
purpose of this section. 

An illustrative example of that situation is the problem of diffusion in a 2D 
channel of non-constant cross-section in the presence of gravity. The situation 
is depicted schematically in Fig. If gravity is involved, the 2D equilibrium 
distribution is given by the usual Boltzmann's factor 

Peq{x,y)^e^ (40) 
Considering the normalization condition 

1 = j Peqix,y)dxdy ^ J Peq{x)dx (41) 

and using the expression (|4o|), one can identify the reduced ID equilibrium 
distribution as 



Peq{x) 



p(-) (~mgy\ kT ( -mgy^\ ( ( 



—2mgw{x) 
kT 
(42) 

where yi(x) and y2{x) define the boundaries of the channel and 2w{x) — y2{x) ~ 
2/1 (x) is its width. The dynamic equation then follows after introducing Peq{x) 



in Eq. (10), yielding 



dP d /„, .dP D(x)mq , , / mqw(x)\ „\ , , 

m = Tx v'^^^^x + 

for the case of a symmetric channel for which y2{x) = —yi{x). 

From the previous expressions one can distinguish different limiting cases: 
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1. In the limit e = "^^Jf^^^ <C 1, we can expand the reduced equilibrium 
distribution in terms of e 

kT I /n/ \\ II 2mg'u;(a;) ^ 2 



^e,(x) ^ —(1 + 0(e)) (^1 - 1 + fcy' ' + O(e^) j ^ 2«;(x), (44) 

recovering the case of purely entropic barrier discussed in the previous 
section. Eq. ( ^3| ) reduces then to the Pick- Jacobs equation (|3C|). 

2. For the case e 3> 1, gravity dominates, cothe ^ 1, and we recover the 
dynamics for diffusion along a ID purely energetic barrier. 

dP d f , ,dP D(x)mg \ 

3. If the boundaries are flat, that is yi{x) = —y2{x) = cte, the equilibrium 
distribution Peq{x) is then constant, and the dynamics can be properly 
described by the ID Smoluchowski equation. 

dP d f , ,dP 



The evolution of this confined system in the presence of gravity dictated by Eq. 
( p3| ) presents some peculiarities. On one hand, since the energy barrier depends 
on the coordinate we have eliminated, there exists a coupling between entropy 
and energy barriers and not a mere superposition of the drifts related to each of 
them, which is reflected by the presence of the term coth ^ '"g^'^(^) ^ ^jjg drift of 
Eq. (p3|). On the other hand, it is important to highlight that Eq. ( ^3| ) does not 
satisfy the detailed balance condition: the mobility b = ^^r- coth ^ ™9'^(^) ^ ^nd 

the diffusion D{x) are not related through the usual Einstein relation D{x) — 
kTb. That is the signature that the fluctuation-dissipation theorem, which holds 
at equilibrium, may lose its validity when we perform a reduction of the variables 
describing the state of the system. 

Concerning the effective diffusion coefflcient, a scaling law as the one pro- 



posed through Eq. (|32| ), which is valid in absence of external forces, will be in 
general not correct when an external fleld is present. The value of the scaling 
exponent a varies in the presence of external potentials, since the scaling ex- 
ponents changes due to interactions. The range of validity of the Fick- Jacobs 
description and the modiflcations in the scaling of the effective diffusion coeffi- 
cient under the presence of external flelds requires a more elaborated treatment. 



6 Conclusions 

In this paper we have presented a theory to describe the kinetics of a system 
whose equilibrium state is characterized by a given landscape of an unspecifled 
thermodynamic potential. The theory is based on mesoscopic non-equilibrium 
thermodynamics, which uses the scheme of nonequilibrium thermodynamics 
at the mesoscopic level of description in which the pertinent flelds are probability 
densities. 
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Particularly, we have established the kinetic equation for a system in the 
presence of entropy barriers. The barriers may be inherent to the intimate 
structure of the system or may emerge as a consequence of the elimination pro- 
cedure of some coordinates when one tries to simplify its dynamic description. 

An illustrative example treated in the literature is the diffusion of a particle 
in an region of irregular geometry. The governing equation for the probability 
density is known as the Fick- Jacobs equation and has been derived heuristically 
by Jacob |P3| and directly from the proper coordinate reduction procedure by 
Zwanzig ||2j]. In our derivation, the entropic barrier directly comes out from 
Boltzmann's entropy, which follows from the proper accounting of the number 
of accessible states of the system. In the kinetic equation we have obtained, 
the dependency of the diffusion coefficient on the coordinate follows from the 
general dependency of the Onsager transport coefficients on the state variable, 
according to the rules of nonequilibrium thermodynamics. This dependency 
becomes crucial in order to the reduced description accurately resembles the 
exact solution. We have proposed a scaling law for the diffusion coefficient 
reaching a very good agreement with the exact solution even in the case when 
the original Fick-Jacobs equation does not provide a good approximation. In 
this way, we have put the validity of a Fick-Jacob-like description in a more 
general context. 

The theory we have proposed can also be applied to cases in which an energy 
barrier is also present. We have obtained the corresponding Fokker-Planck 
equation describing its dynamics. Interesting characteristics of this equation are 
the coupling of the entropy and energy barriers and the apparent violation of the 
Fluctuation Dissipation Theorem, resulting from the elimination of variables. 

The theory we have presented is appHcable to a wide variety of systems of 
different nature. Apart from the diffusion problem in the presence of an en- 
tropic barrier we have discussed in this paper, we can quote protein folding p^ , 
glassy systems ||2^, transport of ions ||2^ and macromolecules j29|-||3^ through 
membranes or channels, motion of polymers subjected to rigid constraints p^ , 
protein binding kinetics jsj, drug release (s^, nucleation js?) or polymer 
crystallization ||3^, to mention just a few examples in which the presence of 
entropic barriers becomes relevant. 

There still remain open questions whose answers goes beyond the scope of 
this paper. The main pitfalls arise from the coordinate reduction procedure. 
The question is how to proceed with the coarsening in more complex situations. 
For instance, it will be desirable to analyze the effect of the asymmetries either 
intrinsic to the nature of the landscape or due to the presence of external forces 
which impose a preferred direction, which occurs in many situations of practical 
interest. Another interesting point will be to analyze the behaviour when the 
slope diverges, giving rise to regions in which the particle may become trapped 
thus breaking the ergodicity of the system. The solutions to these problems 
probably demand a more elaborated analysis of the scaling of the diffusion 
coefficient which takes into account these effects. 

Just proceeding systematically along the lines we have established by means 
of our theory one could arrive at a more complete theoretical framework pro- 
viding the dynamic description of systems exhibiting entropic barriers. 
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